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Abstract 

We analyze, in leading and next to leading order of the BFKL equation, the effects of the 
quantization of the singularities of the j-plane, t-channel partial waves due to the imposition 
of appropriate infrared and ultraviolet boundary conditions. We show that the intercepts, 
uj n of the Regge poles, which contribute significantly to the gluon density in the kinematic 
region measured at HERA and which can be calculated in QCD and in a supersymmetric 
extension of QCD, are substantially modified by Beyond the Standard Model (BSM) effects. 
We also develop a physically motivated heuristic model for the infrared boundary condition 
and apply it to the gluon density. We argue that, using this type of model, the analysis of 
present and future low- a; data could allow one to detect supersymmetry at a high energy 
scale. 
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1 Introduction 



The BFKL equation determines the high energy behaviour of the virtual gluon-gluon scat- 
tering amplitude in Regge limit, in which the cms energy, y/s, is much larger than the 
transverse momenta, k, k', of the gluons. It was derived in the fixed coupling constant case 
by resumming all the Feynman diagrams describing gluon-gluon scattering in the leading or 
next-to-leading order. The scattering amplitude displays a scale invariance such that it can 
be described solely by functions of ratios of transverse momenta k and k'. 

A common application [IH2] of the BFKL equation is to use it as an evolution equation 
in rapidity, y and transverse momentum k for some large-rapidiity amplitude, namely 



and the coupling d s runs with transverse momentum. This application is often used as it 
lends itself relatively easily to an extension of the DGLAP formalism to very low values of 
Bjorken-x, where the pure DGLAP formalism is known to break down. 

In order to solve this evolution equation, one requires as input the amplitude at some 
rapidity y for all values of transverse momentum t'. One could naively expect that eq. (II. ip 
could be directly solved for large y and t well in the perturbative region because it is known 
that the BFKL kernel K(t,t') is quasilocal in t, i.e. it diminishes when \t — t'\ is large. 
However, as was carefully investigated e.g. in ref. [TJ, the BFKL equation is not only an 
evolution in y but also in virtualities t, which leads to a substantial diffusion into the low 
transverse momenta region where perturbative QCD cannot be valid. Therefore, the authors 
of ref. PP proposed a modification of the BFKL equation by imposing a low and high cutoff 
in t on the BFKL integral, a procedure which is today widely accepted. However, such a 
cutoff implies that the amplitude actually vanishes below a certain transverse momentum, 
rather than becoming non-perturbative. In fact, due to the growth of the coupling constant 
at small kx, the amplitude could grow in this region and the vanishing of it at the infrared 
cutoff looks unnatural. 

In ref. [I] we proposed to solve eq. (11. ip by the Green function method which does not 
require any cutoff on the BFKL integral. Instead, we assumed that the non-perturbative 
infrared region of QCD imposes a certain phase on the oscillatory parts of the eigenfunc- 
tions at some small transverse momentum. This treatment of the infrared boundary leads 
to a discrete set of eigenvalues, u n , of the BFKL kernel, since only certain values permit 
the construction of eigenfunctions which simultaneously obey these phase conditions at low 
transverse momentum and the large transverse momentum boundary conditions imposed by 
the asymptotic freedom. This is in contrast to the "usual" treatment [IH3] , in which a lower 
transverse momentum cutoff is imposed on the amplitude, i.e. the amplitude is assumed 
to vanish below the cutoff. In our approach, the amplitudes are particularly sensitive to 
the exact values of the discrete u n which are related to the non-perturbative phases, r] n , 





where 
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at the cutoff. These phases are determined, in turn, by the gluon-gluon interactions of the 
non-pertubative QCD, which lead to rich structures below this cutoff. 

A rigid cutoff (either UV or IR) destroys the scale invariance of the BFKL kernel - and 
hence the validity of scale covariant solutions. As pointed out in p], the imposition of such 
cutoffs has no effect on the position of the leading singularity (the exponent of x in structure 
functions) but does affect the pre-factor, which is controlled by the form of this singularity. 
We find that the subleading singularities are also essential in order to obtain a good fit 
to HERA data. In our application of BFKL dynamics, the scale and conformal invariance 
(which is central to the BFKL formalism for fixed coupling) is broken in a contolled way, 
namely only through the running of the coupling and we assume "quasi-conformal" solutions 
in which the exponent of the transverse momentum varies slowly in order to compensate for 
the change in the coupling in an accordance with a generalized DGLAP dynamics.. 

In our previous paper [I] we have shown that HERA Fi data, at low x, can be described 
very well by the gluon density constructed from the discrete spectrum of eigenfunctions of 
the BFKL kernel. The spectrum contained many eigenfunctions, (9(100), with eigenvalues 
u n varying from U\ ~ 0.25 to u n ~ 0.5/n for large n. This first successful confrontation of 
the BFKL formalism [B] with data led to the unexpected question as to whether the HERA 
data are sensitive to the Beyond Standard Model (BSM) effects. These effects, although 
only present at scales that are much higher than the region of HERA data, can nevertheless 
affect the quality of the fits to data since BSM effects change the running of the coupling 
and consequently also substantially change the values of u n . 

This seems somewhat counter-intuitive. One may ask how it can be possible that a fit to 
data at relatively low energies can be sensitive to corrections due to loops of particles whose 
masses are far in excess of those energies. The crucial point is that the above-mentioned 
large transverse momentum boundary conditions, imposed by the running of the coupling, 
can occur at very high scales 0- For tu <~ 0.1 this is already above the scale at which 
one might expect BSM physics to occur. The value of the discrete eigenvalues arise from an 
interplay between these ultraviolet boundary conditions and the infrared boundary condition 
arising from the imposition of an infrared phase. It is in this sense that in our treatment 
of the BFKL formalism there is communication between high and low energy scales. We 
make the assumption that the allowed eigenvalues can be obtained from the BFKL equation 
supplemented by an infrared phase condition in a process-independent way, i.e. without need- 
ing to impose any infrared or ultraviolet cuts on the integration over transverse momentum. 
Once these eigenvalues are calculated their corresponding eigenfunctions are convoluted with 
the necessary impact factors in order to obtain the required amplitudes. In this way, it turns 
out that even though these amplitudes, in accordance with kinematical constraints, never 
involve diffusion into transverse momenta above the threshold for BSM physics, the rapid- 
ity dependence of such amplitudes is affected by the substantial changes in the eigenvalues 
arising as a result of BSM physics. 

To understand how the running of the coupling constant can have such far-reaching con- 
sequences we derive analytically (in Section 2) the main properties of the discrete pomeron 



Note that this means that formally we determine the eigenvalues at asymptotically large initial energies. 
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solution using the LO BFKL equation. This derivation provides a qualitative physical expla- 
nation of the mechanism by which the BSM effects modify the discrete pomeron structures 
and lead to a genuine change of the eigenvalues and eigenfunctions. It also elucidates the 
role of the infrared phases which define the boundary condition and which can be indirectly 
determined from data. This explanation is then carried over into the NLO evaluation, which 
was used for data analysis and was performed numerically. 

As a popular example of BSM effects we have chosen the N=l supersymmetry and modi- 
fied the /3-function and the kernel of the BFKL equation to include the contributions from the 
superpartners. We then describe the full NLO evaluation of the Discrete BFKL Pomerons 
(DP) with collinear resummation [10J. This allows us to show that the eigenvalues, u n , at 
larger n, have a genuine sensitivity to BSM physics because the support of the corresponding 
eigenfunctions extends to very high virtualities. The values of u n are determined (to large 
extent) by the running of a s and the properties of the BFKL kernel in the high virtuality 
regions, where BSM effects dominate and QCD NLO correction are very small. This is 
also the reason why these eigenvalues are not sensitive to a particular choice of the infrared 
boundary. All these properties are discussed in detail in Section 2. 

In Section 3 we then show that it is possible to construct a physically self-consistent 
infrared boundary condition which determines the properties of the gluon density to be in 
agreement with data. Finally we confront the DP gluon density with the HERA F 2 data and 
show that within our model for the infrared boundary condition we obtain indirect evidence 
of a supersymmetric threshold in a multi TeV range. 

In Section 4 we discuss our results with particular emphasis on the role of the universal 
Green function and the momentum conservation. We also discuss the dependence of our fits 
on the choice of the infrared boundary condition. Section 5 presents a summary. 



The forward amplitude for a diffract ive process with rapidity (or rapidiy gap) y is determined 
by the QCD pomeron and may be written 



where Gu(t, t') is the Mellin transform of a universal (i.e. process independent) Green function 



The process dependence enters only through the impact factors, <3> u , $^ at the top and 
bottom of the gluon ladder which depend on the transverse momenta of the gluons and 
may also depend on other kinematic variables. The integral over u in the inverse Mellin 
transform, (12.2)1 . is performed over a contour parallel to the imaginary axis, to the right of 



2 The Discrete BFKL Pomeron 




(2.1) 



G(t,t',y-y') 




(2.2) 
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all singularities of the Green function. This Green function obeys the equation 

wQ u (t,t')-J dt"K{6t s ,t,t")g w {t",t') = 8{t-t'), (2.3) 

which is solved by determining the set of eigenfunctions of the BFKL kernel, subject to 
certain boundary conditions. If we allow the coupling a s to run with t, t', which in LO 
means replacing it by a/ a s (t)a s (t') then, as we explain below, the UV boundary condition 
- namely that the eigenfunctions decay as t — > oo is automatically implemented. The in- 
frared boundary condition is imposed by requiring that the eigenfunctions have some given 
non-perturbative phases at some low value of t q As was first shown in ref. [12] and we 
explain again below, the combination of the UV boundary condition (which determines the 
oscillation phase at t — t c , where the oscillatory behaviour changes to an exponentially de- 
caying behaviour, compatible with a DGLAP analysis in the double logarithmic limit) and 
the infrared phase condition leads (for positive u) to a discrete set of allowed eigenvalues u n 
with their attendant eigenfunctions, so that the Green function may be written as 

g u (t,t) = yl^lMl + JL f° du' f: ' {t ' )fAt) (2.4) 

^— ' i,i — 1,1 .. Viri I i,i — Cj )' —I— re 



oj — oj n 2ixi J_ 00 oj — oj' + ie 

No further cuts from kinematic constraints on t are imposed in the determination of the 
eigenfuntions and eigenvalues. However, as we discuss in detail in section 4, owing to the 
quasi-local nature of the BFKL kernel /C, this Green function is rapidly attenuated for large 
\t — t'\, which means that when it is inserted into the expression (12. ip in order to obtain the 
unintegrated gluon density, the diffusion into values of t substantially above the region of 
support of the impact factor, $ p (t) is highly suppressed, thereby automatically limiting the 
gluon virtuality to be small compared with the incoming energy. 

The infrared non-perturbative phases, i] n , which determine the values of the discrete 
eigenvalues, ou n , are in general oj dependent, but they must lie within a range of magnitude 
7T, so that the quantum number n represents the number of oscillations of the eigenfunction 
between the scale, to, at which the infrared phase condition is applied and the ultra-violet 
scale, t c , at which the oscillatory behaviour becomes an exponentially decaying one. As 
pointed out in ref. [I], an w-dependent infrared phase condition is necessary in order to 
be able to express an impact factor with support only for small t in terms of the discrete 
eigenfunctions, since the frequency of oscillation of these eigenfunctions at small t is always 
below ~ 0.7. In fact, after imposing this non-perturbative phase, the eigenfunctions form 
an almost complete set of functions in the region of comparatively small k ~ 1 — 10 GeV. 

The existence of this set of discrete eigenfunctions is consistent with the known fact that 
in the Regge regime, the amplitude is determined by a set of Regge poles. The imposition 
of infrared phases does not in any way violate the kinematical constraint, but the ensuing 
discrete spectrum of eigenvalues has a very significant effect on the fitting of the results of 
this modified BFKL formalism to HERA data on the structure functions at low- re. 

The continuum contribution (for negative uj) is not significant for sufficiently small values 
of x at any given t. However, for a given x, as t increases these contributions become more 



These phases should be universal (process independent). 
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significant and are essential in order for this formalism to match the double logarithmic limit 
of the DGLAP approach for sufficiently large t. 

The value of t c , at which the oscillatory behaviour converts into an exponentially decaying 
one, increases linearly with eigenvalue number n. For n ^ 3 this occurs at values of t c above 
the scale at which one may expect to see physics beyond the standard model. The running of 
the coupling is therefore affected by the presence of thresholds for such new physics and this 
in turn affects the positions of the discrete eigenvalues u n . For sufficiently low values of x 
the contribution from all but the first two or three eigenfunctions is negligible. However, we 
find that for values of x which are probed at HERA there is a sizable contribution from these 
higher eigenfunctions. This means that even though the transverse momenta do not diffuse 
into regions of t where the particles of new physics are actually produced, the shifts in the 
positions of the eigenvalues due to new physics affects the x-dependence of the unintegrated 
gluon density. Despite the fact that these effects are small, the high quality of the HERA 
data means that the quality of the fit is significantly affected by the possibility of new physics 
at high energies. 

We now show how this works in detail. 



2.1 LO evaluation 

We begin this section by reviewing the argument of [12] which led to a modification of the 
BFKL formalism which gives rise to discrete poles rather than a cut in the w-plane of the t 
channel partial waves. 

We consider the case of the leading order BFKL equation [6] with running coupling also 
taken to leading order so that (for t > 0) 

_ _ C A a s 1 

a s = = -5-7 2.5 

where 

Po = AC~ A ~ 12 "IS" (2 ' 6) 
The Hermitian BFKL kernel may be written as 



y/a.{t)a a (V)Ko{t,lf), 

where 

r dt'JC (t,t')e M ' = xMe ivt (2.7) 



XoW) = 2*(1) ~ * Q U ~ (2-8) 

Note that the parameter v may be real or imaginary for real eigenvalues Xo{ v )- 
The eigenfunctions, g u (t) of this Hermitian kernel obey the eigenvalue equation 

<tty/a 8 (t)a a (V)Ko(t,l!)g u (l!) = ug u (t) (2.9) 
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These eigenfunctions form a complete orthonormal set 



dtg u {t)g*,{t) = 2n6{u-u/), 

assuming (for the moment) a continuous spectrum for the eigenvalues, u. 
They can be obtained by defining a set of functions f u (t): 



Ut) 



which obey the eigenvalue equation 



a s (t) / dt% {t,t')Ut') = uUt) 



(2.10) 



(2.11) 



(2.12) 



Using eq.f l2.5p we have 



dt% (t,t')Ut') = p utf„{t) 



(2.13) 



Taking the Fourier transform 



Ut) 



and using eq.f l2.7p we have a first-order differential equation 



/3 u 



which has a well known solution 



f u (v) = exp 
In this way we obtain 



(3 u 



Xv{v')dv' 



Xo{v)U») 



r(i/2 + zv) , mi 



(l)u 



V2 



TIOJ 



due 



IV t 



r(i/2 - w) 



T(ll2 + iv)~*, x 



l/03aw) 



(!)«/ 



-i l/(0ou>) 



r(l/2 -iv) 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



where the pre-factor is taken such that the normalization condition f)2.10p is obeyed. 



The integral over v can be performed numerically over a suitable contour. A very good 
approximation to this integral (for small u) is obtained by the saddle-point approximation, 
(equivalent to solving eq.f l2.13p using the semi-classical approximation). The saddle point, 
which is a function of t, z/ s (t), is obtained from the solution to 



Xo{v s (t)) = PoWt. 



(2.18) 
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V plane 



Figure 1: Integration contour (blue line) on the v plane for t > t c . The black dot shows the 
position of the saddle point, v s . 



We consider two regions depending on whether t is greater or less than a critical point, t c , 
given by 

AM C = Xo(0) = 4 In 2. (2.19) 

t >t c : 

In this case there is a single saddle-point on the positive imaginary axis, shown in Fig. HJ 



If we define 7 S by 



1 

7 S = - + %v s 



(2.20) 
(2.21) 



then at the saddle-point, j a , is the solution to 

Xo{ls) = 2*(1) - *( Ta ) - *(1 - 7a ) = o ut 
7 S is in the range 

The contour of integration is deformed so that it becomes the contour of steepest descent 
obtained from the solution to 



arg < / x{v')dv' - x{^s) iy ~ v a ) 



TV 

2 



Near the saddle point the contour runs parallel to the real axis but for very large \u\ it 
runs parallel to the imaginary axis. In the saddle-point approximation, we obtain (using 
eq.(E2TJ)) 

3 ( 7 .-l/2)g(l) r (7s) 



Ut) 



2x'o(% 



r(i-7.) 



V(/9o") 



(2.22) 



V plane 



Figure 2: Integration contour (blue line) on the v plane for t < t c . The black dots show the 
positions of the saddle points, ±|z/ s |. 

This is an exponentially decreasing function of t. Moreover ^ s can be related to the anomalous 
dimension in the DGLAP formalism since 

j t (e t/2 Ut)) = Is (e t/2 Ut)) ■ (2.23) 
From eq. (I2.18P , the anomalous dimension is 

7s ~ | o ( "^) 2 

to \ u 2 

in agreement with DGLAP for small 7 S . 

t < t e : 

Here we have two saddle points lying on the real axis at ±i/ a , shown in Fig. [2j The positions 
of the saddle points are obtained from 

2*(1) - 2Ue |* Q + zi/,J | = (3 ut. (2.24) 

We need to integrate around both of these saddle points, taking a contour of steepest descent 
in the vicinity of the saddle-points, which in this case is inclined at an angle of ±7r/4 to the 
real axis and enclose the positive imaginary axis at large v. The saddle-point approximation 
then yields 

^SK^+f^i)' (2 - 25) 



where 
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The inclusion of 7r/4 in the phase in eg. (12.251) ensures a matching of the solutions at t — t, 
Near t = t c the solution is given by an Airy function. 



€■ 



For t < t c , we have an oscillatory solution which does not lend itself to a match to the 
DGLAP formalism - in this regime a DGLAP analysis is not appropriate, since in this region 
the saddle-point 7 S is complex and double-valued, i.e. 

1 

7, = -±i\v s \, 

and cannot be related to the (real) anomalous dimension of the DGLAP formalism. In this 
region the DGLAP equation is not valid and the BFKL equation can be considered as a 
generalized (quantized) version of the DGLAP equation. 

Recall that the saddle-point, z/ s (t), is a function of t and so we do not have constant 
frequency oscillations. As t — > 0, v s tends to a value vq ~ 0.635 and we have constant 
frequency oscillations in the infrared limit. As t increases this frequency decreases, becoming 
zero at t — t c . 

The infrared phase at t — 0, calculated from perturbative QCD (with t > 0), is then 
given by 

^(^o) . 7T 

Vo* = ~g — + 7, (2-27) 

(0(z/ o ) ~ 0.96). This phase is only determined up to a multiple of tt. We now make a very 
general assumption that the infrared properties of QCD fixes this phase (in general as a 
function of u) to be rj(u), where the function t](uj) is determined from the non-perturbative 
regime of QCD (t < 0). The matching of the two phases 7/0 and 7](u) in the semi-classical 
solution, eq. fl2.25p . then restricts the allowed values of u; to a discrete set u n that satisfy the 
equation 

= vM+(n-\) (n= 1,2,3-..) (2.28) 

The function rj(u) could be a constant (as originally proposed in [12]) but in general it can 
vary with oj. Because of periodicity it can take values in the interval between 0.25 and —0.75 
only. Although i](u) cannot be determined from the perturbative analysis described here, 
its restricted range limits its effect on the determination of the eigenvalues (see Section 2.5). 
However, its variation with u is very important in the construction of the gluon density, (see 
Section 3). 

The above analysis shows clearly that the solution of the BFKL equation has to be 
given by the set of discrete eigenfunctions, whose support in the virtual gluon transverse 
momentum is determined by the critical point, t c , eq. (I2.19p . and whose phase, r](u), at some 
low transverse momentum is determined by the non-perturbative sector of QCD. Effectively 
the boundary conditions modify the BFKL kernel so that it may be written as 

/C(M') = J> n / n *(*)/ n (i') (2.29) 

n 

This coincides with the (LO) kernel 

a s (t)K (t,t') 
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provided it acts on a function, f(t) which may be written as a superposition of the eigen- 
functions f n (t): 

f(t) = J2 a nfn(t)- 

n 

In this way we have supplemented the kernel with both infrared and ultraviolet boundary 
conditions. The infrared boundary conditions arise from the non-perturbative sector of 
QCD, but the ultraviolet boundary conditions arise naturally from the asymptotic freedom 
of QCD. Importantly, the behaviour of the eigenfunctions in the ultraviolet is controlled by a 
critical value, t c , of transverse momentum, which grows almost linearly with n in accordance 
with the fact that the period of oscillations is practically independent on n. Therefore, the 
n-dependent boundary condition leads to qualitatively different results from those obtained 
using a kernel in which boundary conditions are effected simply by a cutoff on a wave 
function. 

Furthermore we note that the value of this critical transverse momentum depends almost 
entirely on the eigenvalue, u, which decreases like 1/n, for large n, as the quantum number 
n increases, (see eq. (I2.28p . This means that in turn the value of the critical transverse 
momenta, k c , increases exponentially as u decreases, so that the n th critical momentum is 
given by (inserting eq.f l2.28p into eq. fl2.19p ) 

e } = A ^ ex p{^y( n -i + ^»))} « A-qcd e 4,5n (2.30) 

Finally let us note that our solution of the BFKL equation is similar to the WKB method 
for the bound-state solution of the Schrodinger equation in the semi-classical approximation; 
the critical point, t c , is analogous to the turning point x c where the potential is equal to 
the energy. Inside a potential well the solutions are oscillatory and outside they decay 
exponentially. This shows that the solution of the BFKL equation (11.11) consists of the 
superposition of the bound state eigenfunctions of the two gluon system with pseudo-energies 
given by the eigenvalues u n . Knowledge of the eigenvalues gives important information about 
the interactions between gluons, both in the infrared and ultraviolet regions of k. We note 
that in the solution of the BFKL equation, the oscillations of the eigenfunctions at large 
k should cancel each other in accordance with the kinematical constraints provided by the 
beam energy of the experiment. This imposes additional restrictions on the non-perturbative 
phases r](u), see below. 

2.2 Threshold effects 

The above analysis assumes that /3o is a constant so that the coupling a s (t) is given simply 
by eq. (12.51) . However, we know that there are thresholds at t — t{ where heavy flavour quarks 
can be produced, and also there may be extra thresholds arising from BSM physics with a 
threshold ( that according to eq. (12.301) can be large) below t = t c . This means that eq. fl2.17l) 
can only be used as a solution for /^(t) between thresholds. As an example, suppose that 
there is only one threshold, at t — t t below the critical point, t c , and that (3 takes the value 
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/?o > above this threshold and (3q < below. At t >t t we have 



a/ 27r/3 w 
and for t < t t we have 

= A J dve ivt 
with the constant A chosen to be 



dv e 



iut 



r(i/2 + ^) 

T(l/2 — iz/)' 



i/(A) > w) 



(2.31) 



T(l/2 + iv) 
r(l/2- zV)' 



-2i*(l)i/ 



i/(A) < «) 



(2.32) 



A- 1 



di/e 



r(i/2 



21/ ) 



r(i/2 - ii/) 



1 !/(/<)") 



(2.33) 



so that the solutions match at t — t f 



In the saddle-point approximation, we can handle such thresholds by noting that eq.f l2.17p 
can be written as 

M) 1 



V2 



ITU) 



(2.34) 



where the "action" S(u, t) is given by 

1 f 

S(v,t)=uvt- -=- / Xo(v')dv'. 
Po Jo 

At the saddle-point, v = v s (t), upon integrating by parts this may be rewritten as 

S(t) = u f u s (t')dt' 



(2.35) 



(2.36) 



where the function u s (t) is given by eq. fl2.18p . (We have used the notation S(t) to denote 
S(u s (t),t) - it is now a function of t only). Here we see explicitly that the saddle-point 
approximation for integral (I2.34p is equivalent to the semi-classical approximation. Replacing 
the integral over v in eq. (I2.34j) by the value of the integrand at the saddle-point we obtain 
a solution which obeys the differential equation 



j t Ut) = x -\MUt) 



(2.37) 



The semi-classical approximation consists of the assumption that the solution to eq. (12.13j) is 
the solution to eq. f!2.37|) multiplied by a slowly varying factor, which turns out to be the same 
as that obtained in the Gaussian integral around the saddle-point in eq. fl2.34p . In analogy 
with the WKB approximation in the Schrodinger equation, there exists a critical point, t c , 
at which the approximate solution changes from an oscillatory function to an exponentially 
decaying one. The assumption of a slowly varying pre-factor breaks down at this point, but 
the solutions either side of the critical point can be matched using a suitable Airy function. 
It is this matching, together with some property of the behaviour at t — that determines 
the allowed eigenvalues. 
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Thresholds are handled in general by replacing eq. fl2.18p by the more general relation 

XoMt)) = (2.38) 
a s (t) 

a s (t) may now be determined using the /3-function with appropriate thresholds. 
The eigenfunctions for t <C t c now take the form 

and the semiclassical quantization condition on the allowed eigenvalues becomes 

S!0i ^(u n )+n-])n (2.40) 



Comparing this with eq. (l2.28p we see that S(0) is independent of tu 0. However, its value 
clearly depends on the positions of the thresholds in /3 , eqs. (12.311) and (I2.32p . Already for 
n > 3 the value of k c given by eq.( I2.30P exceeds many tens of TeV and this means that 
the spectrum is sensitive to any BSM physics. The BSM effect changes the perturbative 
phase r] if the corresponding threshold, k t , is below the critical point, kc ■ The effect of 
such a threshold can be readily estimated; let us assume that for a given eigenvalue, u, the 
threshold t t = ln(kf / Aq C D ) is the largest threshold, below the critical point t c , and that in 
the range 

t t < t < t c , 

flo takes the value (3q > . The exact solution for the LO BFKL with running coupling at t — tt 
is 

f(tt) oc / dve^ v ' u \ (2.41) 



where 

, M() _( it _ 2 JM) +^3,^(1^)} (2.42) 

We evaluate this integral using the saddle-point method and assuming that the threshold 
occurs sufficiently close to the critical point t c that v s is sufficiently small for the diffusion 
approximation to be valid. The phase difference between the case where there is a threshold 
at t = tt and the case where there is no threshold is then 

1 / -< \3/2 1 / X 3/2" 

— (41b 2 - ^utt) - -> 4 In 2 - ^utt 



A0 



3^Vl4C(3) 



o 



(2.43) 



where /3 denotes the value of /3 below the threshold. 



We observe that the above phase difference is substantial for oj < 0.1 and that this 
difference is not suppressed by the scale of the BSM physics. In this sense our analysis 
differs fundamentally from the treatment of the Standard Model as a low-energy effective 
theory way below the thresholds of new physics. In the latter case, logarithmic corrections 
can always be absorbed into unphysical renormalization constants of the renormalizable 
operators of the effective theory, leaving only higher dimension operators whose coefficients 
are suppressed by powers of the new-physics mass scale. 

3 For sufficiently large n we can see this directly since S(0) ~ O^nt^) and whereas tc oc n, uj n oc 1/n. 
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2.3 NLO evaluation 



We have shown in [I] that the BFKL integral equation can be cast in the form of the 
pseudo-differential equation 

d s (t) J dt'Ko(t,t')Ut') = x(-if t ,a s (tyjU(t) = uUt). (2.44) 

Making the simplifying assumption that d s is given by eq. (12. 5p and taking into account the 
collinear resummation jT0j we can write the BFKL equation in the next-to-leading order as 



A)W*/ w (t) 

where the operator v = —id/dt, 
and 



tpo 



2*(i)-*(i+»i/ + |) -* 



Ut), 



1 UJ 

iv H — 

2 2 



(2.45) 



iiy) = X x{v) + - 



2*(1) - * ( - +%v 



* I - - if 



W( ~ + iz> I ( if 



where xi denotes the NLO characteristic function |9J. The function £(z/) does not have the 
poles of the third order of the form ~ l/(l/2±iz/) 3 , in agreement the a renormlaization 
group analysis [9]. 



The equation (I2.45[) can be considered as a quadratic equation in t 

[u(tp ) 2 -tMo(i>, ")-&)] ■ f u (t) = 0. 



(2.46) 



We can convert this into a second order differential equation for the Fourier transform fcj(v) 



where t is replaced by the operator t 



' du ' 



In the semi-classical approximation in which 



ln(xo) and ln(£) are treated as slowly varying functions of v, so that 



2 _ 



this second-order differential equation may be written as 

[iM - X~{v, w)] • [0 o t - X + (u, w)] • f u {y) = 

where 

X ± (v,u) 



(2.47) 



1 



2u 



1 

2oo 



with solution for f u (t) which is analogous to the LO expressions of eq. (12. 14|) and (12. 16^ : 



Ut) 



dve lvt ex$ ( -j- / X + (u,cu)du ) . 



(2.48) 
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For small u, eg. ( 12. 48 ft may be approximated by 



Ut) = -= / cfc/e 4 * 
V27rw 



-oo 



T(l/2 + iv) __ 2mi 
T(l/2 -iu) 



l/(/8ow) 



e -^r-'^')_ (249) 



Note that the saddle-point of the integral over v in eg. (I2.49P occurs at v s where the RHS of 
eg. (12. 45 P vanishes as reguired for f u (v s ) to be a turning-point. 

In a more general NLO approach, the BFKL eguation (12.441) can be simplified using the 
semi-classical approximation, i.e. assuming that the t-dependence of In f u (t) is large so that 

" y Ut)«Ut)(^W. (2.50) 



dt J \ dt 

The eg. (I2.44p looks then like the non-linear differential eguation 



cHnjU^ = x{]/{fhnAf)) = ^ (2.51) 



dt 

As a result the freguency v(t) is a function of t such that 

( a s {t)C A \ ( a s (t)C A \ 2 

u = I — - — JxoO) + ( — - — J xi\y) + --- ( 2 - 52 ) 

The expression ( 12. 52ft -including collinear resummation [10] - is the NLO analogue of the eg. 
(12T38|) . The eg. (l23T]) has a solution 

/«(*) = e iS(t)/ " (2-53) 

where 

5(t) = u I v{t')dt', (2.54) 

The critical logarithmic transverse momentum, t c , is the value of t for which u(t) = 0. This 
condition is the NLO analog of eq.( l2.19j) . For t < t c , there are two real solutions for u(t) 
generating an oscillatory solution with a given phase, whereas for t > t c the solution is on 
the positive imaginary axis, generating an exponentially decaying function as t — > oo 

Thus we see, that the solution for the eigenfunctions in semi-classical approximation 
is analogous to that in leading order, but the function v(t) takes into account the NLO 
characteristic function as well as the NLO running of the coupling and the threshold effects. 
A further feature of threshold effects beyond leading order is that it is not only the /3-function 
that has steps at the thresholds but also the NLO contributions to the characteristic function, 
5xi, - corresponding to the presence of new particles at some point in the ladders 

The semi-classical approximation is valid provided 

dln{u{t)) 



dt 



< u{t) 
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This condition breaks down in the region t ~ t c where \v(t)\ is very small. However, in this 
region the eigenvalue equation approximates to Airy's equation with solution 



/«(*) = Ai U— > ( 2 - 55 ) 




For t 3> t c the Airy function, Ai, behaves as 



AiH^y 7 j ^ e -m)\/" (2.56) 



and for t <C £ c 



(2.57) 



We therefore find that the solution eq.f l2.55p is a good approximation over the entire range 
of t and at the same time determines the phase of the oscillatory solution for t = t c required 
to match the oscillatory region and the exponentially decaying region. As in the LO case, 
we make a very general assumption that the infrared (non-perturbative) properties of QCD 
impose some phase, r], at t — 0, defined up to an ambiguity of nn, which can also depend 
on uj. We find then that we can only match this phase to the perturbative one, determined 
from eq.f l2.54p . for one value of u for each integer n, where n corresponds to the number of 
oscillations. This leads to the quantization of the spectrum (i.e. discrete pomeron poles) 
given by eq. fl2.40p . in keeping with the predictions of Regge theory. 

In contrast to the LO evaluation, in full NLO the eigenvalues and eigenfunctions can only 
be determined using numerical methods of iteration and integration. Their construction 
requires several steps. In the first step we determine the values of the frequency v as a 
function of u and t from the solutions of eq.f l2.52p . Then, the critical point, t c , is determined 
as a function of u from the condition u{t c ) = 0. The phase function S(t), for a given oj, is 
then found from eq. (I2.54p . 

In the next step the phase, r], at the infrared boundary has to be specified. In the 
leading order computation it was possible to define it at Aqcd, because the frequency v is 
well defined at t — 0, eq.f l2.18p . For the NLO calculation, we obtain v with the help of 
eq. fl2.52p . which is not valid at Aqcd- We therefore defined it as a phase condition at the 
lowest possible value of the (logarithmic) transverse momentum, t = i , which can be safely 
reached by the perturbative calculation (see also the discussion in Section 3.2). 



2.4 N=l Super symmetry at Various Thresholds 

We have chosen as example of "new physics" the popular N = 1 supersymmetric extension 
of the Standard model above a given threshold in kx, which for simplicity we assume to be 
a common mass threshold for all superpartners. Below this threshold the running of the 
coupling is governed by the /3-function to two-loop order 

a 2 s (IIC A 2 \ a 3 s (MC\ (10C A ^ \ \ fn C£> , 

fj < = -47(^-H^(V + (^ +2 ^J % J' (2 - 58) 
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Figure 3: The running of a s across a threshold for N=l SUSY at 10 TeV 

where for the case of QCD, Ca = 3, Cp = 4/3 and rif is the number of active flavours. 
Above the threshold, the /3-function is given by 

P> = "if (3Ca - %) - j^y 2 (eCi + (-^ + 2C F ) n/ ) • (2.59) 

This leads to a "kink" (discontinuity in the derivative) in the running of a s at the threshold 
for N=l SUSY, which can be seen in FigJSJ 

Furthermore, above the SUSY threshold, the NLO characteristic function, Xi{ v ) acquires 
an additional contribution El of 



, x 7T 2 smhi-Kv) ( 11 \ 



from the octet of Majorana fermions (gluinos), and 

r t \ 71-2 n f sinh(7rz/) /5 2 \ , n 

from the squarks (in the fundamental representation). 

Typical graphs contributing to 5 s Xi{v) are shown in Fig. HJ They contribute only at 
NLO level since the exchange of a fermion or scalar particle in the t-channel is suppressed 
in LLA [9], and therefore only contributes at subleading logarithm order. 



2.5 The Discrete Pomeron with and without SUSY 

In this sub-section we investigate the properties of the Discrete BFKL Pomeron with and 
without SUSY contributions. For this example, we have assumed that the SUSY threshold is 
at 10 TeV. Fig. [S] shows the spectrum of the eigenvalues u n computed in the NLO computa- 
tion assuming that the non-perturbative phase, r)(u n ), of the eigenvalue condition, eq. fl2.40p . 
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(a) (b) 
Figure 4: Typical graphs contributing to BFKL kernel involving (a) gluinos or (b) squarks 
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Figure 5: The eigenvalues computed in the NLO evaluation of the Standard Model (triangles) 
and SUSY at a threshold of 10 TeV (squares). The lines indicate the maximal possible spread 
due to the uncertainty of the phase (77) choice. 

is 77 = —0.25, for all eigenfunctions. The perturbative phase, 5(0), of eq.( l2.40p . at the infrared 
boundary t — 0, is replaced by S(k ), with k = Aqcd exp(t /2) and k = 0.6 GeV. The 
eigenvalues determined with and without SUSY effects differ substantially for n > 3 whereas 
for n < 3 they show no difference. This is understandable from the Appelquist-Carrazone 
theorem [7] because the assumed SUSY threshold, that we have chosen, lies between the 
critical momenta for the second eigenfunction {k c ~ 1 TeV) and the third eigenfunction 
{k c ~ 100 TeV). The k c values computed at NLO, with and without the SUSY threshold, 
are shown in Fig. [HI They turn out to be very close to the leading order values calculated 
from eq. (12.301) - the difference being due to the fact that a s runs more rapidly for NLO than 
for LO. Furthermore these critical momenta show only small dependence on the presence on 
SUSY threshold. 

On the other hand, the eigenvalues, u n , which are important for the description of the 
HERA structure-function data at low-x, are very sensitive to possible threshold effects, i.e. 
they differ substantially already in LO and the differences are much larger than any possible 
uncertainties due to the unknown phase 77. For example, in LO, eq.f l2.28p . the ratio of /Vs 
below and above the SUSY threshold is 7/3, which means that already for n > 3 the effect 
of the change in /3 on the eigenvalues, is much larger than the maximal possible effects due 
to the uncertainty in 77 (rj can only vary between rj = 0.25 and 77 = —0.75). At NLO, we 
find that the phase independent discrepancy between the eigenvalues with and without the 
SUSY threshold starts at n > 5 (see Fig, ED . These substantial differences are related to the 
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Figure 6: The critical point k c computed in the NLO evaluation of the Standard Model 
(triangles) and SUSY at a threshold of 10 TeV (squares). The red lines show the LO 
computation which is not sensitive to any threshold effects. 
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Figure 7: The eigenvalues computed in the NLO evaluation of the Standard Model (triangles) 
and SUSY at a threshold of 10 TeV (squares). The lines show the LO computation in the 
two cases. 

fact that for u < 0.1, the change of phase of an eigenvalue, A0, arising from the change in 
Po as one crosses the SUSY threshold is large (as can be seen from eq. fl2.43p ) and indeed 
much larger than the maximal possible rj change, Arj < 1. 

Since the properties of the eigenvalues are determined by the behaviour at very high 
virtualities, (of the order of k c ), it should be expected that the eigenvalues computed in NLO 
should approach the LO ones at large n Q . Fig. [7J shows the comparison of the eigenvalues 
computed using the NLO and LO approximations and confirms this expectation. The LO 
computation was made using eq. fl2.28p with f3 values computed with rif = 6 below the 
SUSY threshold of 10 TeV. It is interesting to observe that LO and NLO results approach 
each other more slowly in the case of the SM+SUSY than in the SM alone. This slower 
approach is due to the fact that a s runs more slowly above the SUSY threshold. This means 
that the eigenvalues, u n , approach zero at a different pace, as can be seen from the figure. 
We note that for small u, the eigenvalues are very closely packed and so the effect of the 
discrete nature of the solutions becomes less important. To a good approximation we could 
replace the sum over the eigenf unctions for small u by an integral over a range of small uj. 
However, it is important to note that the jacobian for the transition from a discrete sum to 
an integral is proportional to the gradient of the u — n distribution shown in Fig. [7J and this 
is different in the two cases - leading to different pomeron amplitudes. 



4 For sufficiently small uj, the NLO effects both in /?o and the characteristic function, %, become negligible. 
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Figure 8: Oscillation frequencies as a function of gluon transverse momentum for various 
eigenfunctions. The left-hand panel is the case of the Standard Model and the right-hand 
panel is the case of N=l SUSY above a threshold of 10 TeV. For the purpose of this com- 
parison it has been assumed that the infrared phases are the same in both cases. 



In addition to the change of the running of a s , there are also effects due to the 5\i 
contributions to xi which sets in above the SUSY threshold, eq. (12.601) and eq. fl2.6ip . It is this 
discontinuity which is responsible for the discontinuities in the frequencies v at threshold, see 
Fig. and not the change in the rate of running of the coupling, which remains a continuous 
function 0. The change in frequency thus compensates for the change in the characteristic 
function in order to ensure that the eigenvalues, u n , remain unchanged as one passes through 
the threshold [§. The contribution, Sxi, of these additional terms is shown as a function of 
frequency in Fig. [9] where it can be seen that this is a rapidly decreasing function, which 
explains why the discontinuities in frequency at threshold are much larger for the lower 
eigenfunctions for which the frequency at threshold is lower. 

For lower n eigenfunctions, the change of the frequencies due to the SUSY threshold leads 
also to the change of its shape. In Fig. [10] we show a representative subset of eigenfunctions 
in the Standard Model and the SUSY model in the transverse momentum region relevant for 
a fit to HERA data. The eigenfunctions are shown with values of rj = —0.25 and rj = +0.25 
(in order to give an indication of the sensitivity of the allowed eigenvalues to the unknown 
infrared phases), with and without SUSY at a threshold of 10 TeV. As expected, the first 
two eigenfunctions are identical since their values of k c are below the SUSY threshold. The 
third and higher eigenfunctions display significant differences for both displayed 77 values. 
Remarkably, these differences diminish for higher eigenfunctions and for n > 50 the two 
eigenfunctions almost overlap in the displayed kx region (relevant for a fit to HERA data). 
The reason for this can be seen from Fig. [HJ which shows that for the relatively low transverse 

5 A similar smaller discontinuity can be seen at around 3 GeV. This corresponds to the c-quark threshold. 
There are analogous, even smaller, discontinuities at the b-quark and t-quarks thresholds 

6 The discontinuous changes in frequency are due to the fact that the change in characteristic function is 
imposed at a threshold in its entirety. A determination of the NLO characteristic function which accounted 
for the mass of internal particles would smooth out these discontinuities. 
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Figure 9: Decrease in the NLO characteristic function, xi as a function of frequency u 

momenta the differences in the frequencies between the two models decreases with increasing 
eigenvalue number, so that if the infrared phases are equal, the functions will be almost 
identical in this region. 

In summary we can state that the Discrete BFKL Pomeron shows a clear sensitivity to 
BSM physics effects and that these effects cannot be absorbed into its only free parameters, 
the infrared phase i](u). This is clearly seen in the eigenvalue dependence on the SUSY 
thresholds, both in the LO analytical approach, eq.f l2.28p and eq.f l2.43p . and in the NLO 
numerical evaluation, Fig. |5j In addition, the asymptotic behaviour of the eigenvalues with 
increasing n (or decreasing oS) is very different for SM and SM+SUSY irrespective of the 
possible higher order QCD corrections, Fig. This means that we have here a very different 
situation from the scenario described by the decoupling theorem [7] where the large loga- 
rithmic corrections can be absorbed into unphysical renormalization constants leaving only 
higher dimension operators whose coefficients are suppressed by powers of the new-physics 
mass scale. In the case of the DP the effects of SUSY thresholds produce large changes of 
frequencies, Fig. [HJ which modify the infrared perturbative phases rj. These alter, in turn, 
the spectrum u n and hence the properties of the gluon density. Since the gluon density is 
a measurable quantity the non-perturbative phases, 77, can also be measured although only 
indirectly. 



3 Comparison with HERA Data 
3.1 General considerations 

One of the most important results of the HERA experiments is the measurement of the gluon 
density. This density encompasses the properties of the pomeron in the sense that the same 
gluon density determines the dynamics of the inclusive 7*p (or F2) and diffractive processes, 
in particular the exclusive vector meson production. Several investigations performed in the 
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Figure 10: Comparison of a representative subset of eigenfunctions in the Standard Model 
(blue) and the SUSY model (red) computed at rj — —0.25 (solid line) and r\ = +0.25 (dashed 
line). The SUSY threshold is assumed at 10 TeV. The eigenvalue number is given in the 
upper right corner. 
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context of the dipole models [T8H25] have shown that the effective intercept of the gluon 
density measured by the rise of F2 with diminishing x, called A, is properly translated by the 
optical theorem, to the effective intercepts seen in the exclusive vector meson production. 

The effective intercept A measured at HERA varies between A ~ 0.2 at Q 2 = 10 GeV 2 
to about A = 0.35 at Q 2 = 100 GeV 2 , see Fig. 9 of [4]. The Q 2 dependence of A in F 2 and 
in its diffractive counterparts can be well reproduced by the DGLAP evolution in which the 
values of A are almost entirely of perturbative origin. In the well-known Donnachie-Landshoff 
(DL) [27] picture of the Pomeron the variation of A with Q 2 is due to the existence of a hard 
(A = 0.4) and a soft, non-perturbative, (A = 0.08) Pomeron which ad-mixtures varies with 
Q 2 - 

The properties of the gluon density corresponding to the DP are determined by the Green 
function constructed from the discrete eigenfunctions of the BFKL kernel (convoluted with 
the proton impact factor). In contrast to the DL Pomeron the DP is composed of infinite 
many eigenfunctions with eigenvalues varying like u n ~ 0.5/n. The eigenvalues u n are almost 
entirely of perturbative origin because its only non-perturbative ingredients are the infrared 
phases rj n , which have a negligible importance for larger n, as was explained in the previous 
section. 

The infrared phases have, however, a strong influence on the shape of the gluon density 
since they determine how the contributing eigenfunctions add together. Let us recall that 
the un-integrated gluon density from the DP is of the form 

g(x, k 2 ) DP =-k 2 T(-) A n f n (k, tin) (3.1) 
x ^— \ \x I 

n=l x ' 

where g means differentiation of the gluon density w.r.t. ln(A; 2 ). Here, the eigenfunctions as 
a function of k (rather than t) are normalized w.r.t. k and are related to f n (t) by 

fn(k,Vn) = ^^Aj^))' (3 ' 2) 

A n is the overlap integral of these eigenfunctions with the proton impact factor, eq. fl3.3p . 
with the eigenfunctions f n (k, rj n ) computed with a specific r\ — n relation, r\ n (after accounting 
for the non-zero overlaps of the eigenfunctions of the non-Hermitian kernel, - for details 
see |1]). The sum over eigenfunctions in eq. (13. 1 p is limited for numerical reasons (see below) 
to n mnT = 0(100)0. The oscillation frequencies of the eigenfunctions, at transverse momenta 
relevant to HERA, vary very little from one eigenfunction to the next, so in order to obtain 
a positive gluon density, which grows with k, it is necessary to generate a strong 77 — n (or 
equivalently rj—u) dependence; the eigenfunctions oscillate in In k and the only way to cancel 
these oscillations is to introduce a shift of the phase between the different eigenfunctions. 

7 Wc have shown in ref. [¥| that an adequate description of HERA F 2 data requires O(100) eigenfunctions. 
Nevertheless, the limit of n max ~ 100 represents a model assumption which is sufficient for HERA data but 
which could be too low for evaluation of LHC DY data. For LHC data it could also be necessary to include 
the contributions of negative w's, eq (|2.4|) . 
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The sum of eq. (13. ip determines the evolution properties of the gluon density in agreement 
with the BFKL equation (11.11) . The infrared phases are determined by the dynamics of 
non-perturbative QCD and it should be possible, in principle, to estimate them using lattice 
calculations. It should also be possible to determine them directly from data at comparatively 
small Q 2 . However, in order to be able to extract these phases accurately, the required 
data set should have a much larger x and Q 2 range than the presently available HERA 
measurements. In addition the data set has to achieve the precision of the present HERA 
F 2 data. 

Therefore at present, to be able to confront the DP with data, we have no other choice 
but to construct a heuristic model for the infrared boundary condition based on the general 
understanding of the non-perturbative physics. We first postulate that the form of the in- 
frared phase function, t)(uj), dictated by the infrared behaviour of QCD, is not sensitive to 
BSM physics and is a smooth function of u rather than an arbitrary number for each eigen- 
function. We then choose to describe it in terms of a suitable parameterization, motivated 
by a similarity of the BFKL dynamics with the Schrodinger equation, described in detail 
below. In addition to the phases, we have also to specify the proton impact factor. 

We consider this heuristic approach as a first step towards the determination of the 
infrared phases, which are important quantities of the non-perturbative QCD. The main 
purpose of the investigation reported in ref. [I] was to check whether a physically plausible 
boundary condition provides a good description of data, i.e whether the Discrete BFKL 
Pomeron can describe the dynamics of the measured gluon density. The main purpose of 
the present investigation is to find out whether the genuine sensitivity of the DP to the 
presence of BSM physics at high-energy can improve or worsen the quality of the fit to data, 
notwithstanding the uncertainties associated with the infrared phase conditions. 

3.2 The Infrared Boundary 

Our heuristic model of the infrared boundary consists of a set of physically well motivated 
assumptions about the proton impact factor and about the rj — n (or equivalently rj — u>) 
relation. 

The proton impact factor has to be positive everywhere and concentrated at the values 
of k < 0(1) GeV. We therefore choose a very simple possible form 

$ p (Jfe) = Ak 2 e- bk \ (3.3) 

as in ref. [I]. We have also investigated other forms of the proton impact factor, e.g. with 
different powers of k 2 in the prefactor and/or the exponent but found that the fit to data 
has no sensitivity to such alternatives. This is due to the fact that all eigenfunctions have 
a similar, oscillatory, shape near the infrared boundary and that the period of oscillations 
of the eigenfunctions is much larger than any physically possible support of the proton 
impact factor. Note that the form (I3.3P vanishes as k 2 for small k, as required by colour 
transparency and that the coefficient b has the interpretation of the average inverse square 
transverse momentum of partons inside the proton (the value of the parameter b was left 
though completely free in the fit). 
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Our choice of ansatz for the dependence of the infrared phases, r/ n , on the eigenfunction 
number, n, is motivated from an examination of eq. (12. 28ft for the eigenvalues at L.O. We see 
that for large n we have u n oc 1/n - the eigenvalues decrease and become closely packed 
as n increases. This is similar to the eigenvalues of a bound state in a Coulomb potential 
problem. The value of r] n has a restricted range (in order to avoid "cross-over" between 
adjacent eigenvalues) and its variation with n must be smaller than n. Since they are 
generated by the quasi-bound states of gluons inside the proton they should be described by 
a simple parameterization. In ref. [3], we found a simple, two-parameter, form 



where n max is the number of eigenf unctions we use for the fit and r] r represents the total 
range (in units of ir). The parameter k must be less than one. (Note that for n max — > oo 
and fixed n, the phase i] n formally tends to 0.) 

Eq.f l3.4n is by no means unique and we could have added terms which are analytic in u n 
of the form 



We have tested such more general parameterizations and found that, despite the introduction 
of extra parameters, there is no improvement in the quality of the fit obtained. We therefore 
use the simple ansatz (13.41) . but we treat rj r as a free parameter (with the only restriction 
that it must not exceed one), in order to assure a bias- free evaluation in all of the fits that 
we perform. 

In ref. jl], we defined the infrared boundary as a phase condition at the lowest possible 
value of the transverse momentum, k = k , which can be safely reached by the perturbative 
calculation. To make this value as close as possible to Aq C d we considered only the one-loop 
running of the coupling. This gave a value of k = 0.3 GeV, which corresponds to a s ~ 0.7. 
The reason for running the coupling at one-loop only was that in principle this is the same 
order of perturbation theory as the NLO characteristic function, xi [H]- However, given that 
we modify the eigenvalue eq. (I2.52H by resumming all the large corrections in xi using the 
technique of ref. [10J . it is more appropriate to take the /3-function to two-loop order which 
is what we use in this paper. 

When we do this, we are faced with a problem - namely that we cannot run the coupling 
below an "infrared" scale ko = 0.6 GeV, which corresponds to a s ~ 0.7 (at the two loop level), 
without approaching the Landau pole too closely. On the other hand, the infrared boundary 
conditions are to be imposed at a transverse momentum of order Aqcd- Moreover, we need 
to know the eigenfunctions below ko in order to perform a convolution with the proton 
impact factor, which has support mainly below ko- Therefore, guided by the behaviour of 
the eigenfunctions in the perturbative region, we continue them down to a lower momentum 
k Q , which should be of order Aq C d, using the extrapolation of the phase (j) n {k) 



\ko/ 

where for each eigenfunction, with index n, z/° is the frequency of the oscillations near 
k = ko [3]. We have assumed that this frequency is constant below k , an assumption which 




(3.4) 



b + cu n + duol H . 




(3.5) 
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is correct for sufficiently small ko, at least for the leading order BFKL kernel (see [12] )■ Any 
deviation from constant frequency should have a negligible effect as we are only extrapolating 
over a small range in gluon transverse momentum. The numerical values of z/° are obtained 
by inverting the eigenvalue equation f)2.52p . modified according to [TO] . 

The overlap integrals between the proton impact factor and the eigenf unctions must also 
start at k (the support of these impact factors being significantly attenuated at ko). We 
therefore use this momentum at which we impose the infrared phases of the eigenfunctions. 
The relation between the phases at fco an d fco is given by eq. fl3.5p . We leave the exact value 
of ko as a free parameter with the restriction that it must be O(Aqcd) and define it to be the 
scale at which the phase of the leading eigenfunction vanishes (as can be seen from eq.f l3.4p ). 

3.3 Results of the fit 

Before a comparison can be made with the measured structure function, F2, it is necessary 
to convolute the gluon density with the impact factor for the virtual photon (for details of 
the procedure see Section 6 of ref. [1]). The impact factor for the virtual photon is calculable 
in perturbative QCD and has support which is peaked at transverse momenta of the order 
of the photon virtuality, \[Q 2 . 

The fits were performed using the HERA data [5] in the low- a; region, x < 0.01. To avoid 
any saturation effects we have limited the fit to the Q 2 > 8 GeV 2 region. We recall that the 
saturation scale at HERA was determined to be Q 2 S = 0.5 GeV 2 at x ~ 10 -3 [2"3"||2"4] . therefore 
our choice of the Q 2 region is very conservative. This choice means that, in this paper, we 
concentrate on the one-pomeron exchange, without any multi-pomeron contributions, which 
could induce saturation effects. The saturation effects could also play a role without multi- 
pomeron effects through a modification of the boundary conditions, see ref. [26], which in 
turn could modify our ansatz for the infrared boundary. In any case our choice of the Q 2 
region for fits assures that saturation effects can be ignored in this first evaluation. In future 
we plan to extend our analysis into the Q 2 region which could be more sensitive to saturation. 

In the region of Q 2 > 8 GeV 2 we have a total of 108 data points and a total of 5 parameters 
- so the number of degrees of freedom is N^f = 103. We consider the Q 2 > 8 GeV 2 region 
as our main investigation region and use the Q 2 > 4 GeV 2 as a cross check. 

As discussed in the previous paper [I], in order to obtain the most accurate estimate of the 
un-integrated gluon density we should include in the fit as many of the higher n eigenfunctions 
as possible. Indeed, we observe that the fit quality improves with increasing number of 
included eigenfunctions and the series converges in x 2 - 111 principle this convergence should 
improve as n — > 00. However, in practice the number of eigenfunctions used in a fit is limited 
by the numerical precision of our calculation. We have indications that if we take significantly 
more than 100 eigenfunctions our fit could be polluted by numerical instabilities arising from 
an accumulation of computational rounding errors. Moreover, we find no improvement in 
the quality of our fits, either in the case of the Standard Model or for MSSM SUSY at any 
of the thresholds investigated, when the maximum number of eigenfunctions, n max , exceeds 
100. We have therefore taken n max = 100 throughout. 
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Table 1: Fits for N=l SUSY at different scales. The bottom row corresponds to the Standard 
Model. All fits are performed with n max = 100. 



In Table CD we show our fits for various SUSY thresholds as well as the Standard Model. 
Let us first note that the ko values obtained in the unbiased fit, k ~ 275 MeV, are close to 
Aqcd- At the same time the value of b implies that the proton impact factor peaks around 
Aqcd, as expected for a self consistent description. This together with the relatively low 
X 2 's of all fits confirms the success of our construction of the infrared boundary. 

The quality of the fits shows a clear preference of the evaluation with SUSY effects; the 
fit for the Standard model is worse than the fits with SUSY thresholds larger than 3 TeV. 
A SUSY threshold of 3 TeV, which is close to the reach of LHC also gives a worse fit. On 
the other hand for a SUSY threshold in the region of 10 - 15 TeV, the quality of the fit is 
significantly improved, but that for significantly larger SUSY thresholds worsens again. 

Let us note that the differences between the \ 2 fits shown in Table [1] are very significant 
because the Maximum-Likelihood method, which assures that the minimum of \ 2 provides 
the best estimate of the parameter values, also states that lcr error in determination of 
parameter values is given by A% 2 = 1, irrespective of the number of degrees of freedom. 
Therefore, the differences between the fits of Table [Hare a multi-cr effect, within our model 
of the infrared boundary. 

The Ax 2 = 1 rule is valid estimate of the parameter error only for estimates 

within one theoretical framework, i.e. for one likelihood function. On the other hand, when 
X 2 is used to quantify the agreement of different theories with data, it is expected that the 
observed x 2 can deviate from the optimal value because it fluctuates with the probability 
density function, f(x 2 , Njf), which is approximately a Gaussian with the average value equal 
to Ndf and the variance, o 2 = 2A^. In the case of fits, presented in the Table 1, the expected 
X 2 should be around 103 and a = 14. Therefore the x 2 values obtained for best fits with 
the SUSY mass 0(10) GeV lie well within one standard deviation. The DGLAP fits have 
a x 2 /Ndf ~ 0.95. This would lead for our sample to x 2 = 98, which is also within one 
standard deviation of the optimal value, so that one cannot conclude that either fit is better. 
The evaluation with the goodness of fit criterion called p-value, which is more appropriate 
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Figure 11: Comparison of the DP fit with M SUSY = 10 TeV with HERA data. 



for Ndf ~ 100, gives a p-value ~ 30% for SUSY masses of 0(10) GeV, which is again an 
excellent result, see ref. [TBI ITT] . 

As a check, we also performed the fits with a lower Q 2 cut, Q 2 > 4 GeV 2 . We find that 
in this Q 2 region there is an significant increase of x 2 /Ndf presumably due to various higher 
order effects, such as the NLO contribution to the photon impact factor or valence quarks 
effects and possibly also the proximity of the saturation region. Although the overall quality 
of the fit for all data with Q 2 > 4 GeV 2 is significantly worse than with Q 2 > 8 GeV 2 the 
preference for N=l SUSY with the threshold region of 10-15 TeV is also clearly seen. In the 
Q 2 > 4 GeV 2 region there are 128 points and the x 2 's of the best fits are 184.3 (3TeV), 
164.5 (6TeV), 155.6 (lOTeV), 152.6 (15TeV), 169.7 (30TeV), 164.7 (50TeV), 164.3 (90TeV). 
The best x 2 f° r Standard Model is 169.7. The values of the fit parameters are similar to the 
values shown in Table HJ for the Q 2 > 8 GeV 2 region. 
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Figure 12: Evolution of the initial wave-packet in the DP solution of the BFKL equation 
as determined from the fit to HERA F 2 data. 



4 Discussion 



The main result of this paper - namely a possible sensitivity of HERA data to BSM effects is 
especially astonishing as it was not seen in the usual solutions to the BFKL equation [TH3] . 
This sensitivity is due to the fact that we use the Green function method, in which the 
Green function is considered to be a universal property of the BFKL equation and that 
apart from supplementing this equation with a given set of infrared boundary conditions 
for the eigenf unctions, no cuts on transverse momentum are applied, so that the conformal 
invariance is broken in a very smooth way, solely by the running of the coupling constant. 

Had we imposed kinematic limits on the gluon virtuality for the individual eigenfunctions, 
then for those eigenfunctions with u J$ 0.1, they would never reach the region in which they 
decay as required for compatibility with a DGLAP analysis in the DLL limit. Alternatively, 
one might attempt to fit HERA data only with the first eigenfuction, which does indeed start 
to decay within the HERA region, making the assumption that the higher eigenfunctions 
(with lower u>) only couply weakly to the proton. However, the properties of this first 
eigenfunction are in clear contradiction with data; for example, the rate of rise of F 2 with 
diminishing x would be independent of Q 2 with a value of A ~ 0.25 whereas the A value 
in data varies between 0.15 to 0.35 in the observed Q 2 region. The results of ref. [UJ 
indicate that one needs at least four eigenfunctions to describe HERA data (making an 
artificial assumption that the overlap constants A n of eq. ( 13. 1 p are totally unconstrained), 
which extends the virtuality region to kf ~ 1000 TeV. 
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In the solution which we have adopted the momentum conservation is incorporated by 
convoluting the Green function with the proton form factor as defined by the solution of the 
BFKLeq. ffTTj) 

/r 
du J — e "yg u (k,k')%(k'), (4.1) 

where in analogy to eq.f l3.ip the Green function, Q w (k, k'), is written in terms of the eigen- 
functions f n (k,r] n ) and f n (k',r] n ). This expression describes the evolution of a wave-packet, 
WV, from rapidity y = to the larger rapidity y values 0. The (approximate) momentum 
conservation emerges here because the quasi-local nature of the kernel )C(t,t') ensures that 
there is no evolution into the very large transverse momenta, apart from the usual BFKL 
diffusion, In kx ~ V a In s. Fig. [12] shows the wave-packet as it evolves from the smallest to 
the largest rapidity values of the HERA region, y = ln(l/x), and for the SM and SM+SUSY 
cases. 

The initial wave-packet, the curve WV(y = 0, k 2 ) shown in Fig. [121 is n °t quite the same 
as the function <& p {k) defined in eq. (I3.3[) (divided by k 2 ) but rather it is related to it by 
eq. (14.11) . i.e. it is built only out of the eigenfunctions which obey the imposed boundary 
conditions on the Green function. For large y, it is sufficient to consider only the discrete 
eigenvalue part of the Green function in eq. (l4.1l) . and so we are actually using that part of 
the proton impact factor which is orthogonal to the (continuum) negative u eigenfunctions. 
The initial wave-packet is somewhat broader than the distribution of eq. (l3.3p but it has the 
required feature that it is localized to a transverse momenta of ~ 1 GeV. With increasing 
rapidity it broadens due to the BFKL diffusion but this broadening remains moderate, of a 
size of just few GeV, as in the usual solutions of the BFKL equation. Although our kernel 
is constructed from the oscillating eigenfunctions with properties determined by the BFKL 
dynamics at very high transverse momenta, the oscillations cancel away due to the choice of 
the phases r\ n in the fit procedure and what remains is only a slight broadening of the gluon 
diffusion spectrum with increasing y. 

In the application to DIS it is more customary to consider the evolution of the uninte- 
grated gluon density given by 

xg(x,k 2 ) = k 2 J du j ^ l^j G u {k,k')%{k'), (4.2) 

which differs from the evolution of the wave-packet, eq.( I4.ip . by the kinematical factor k 2 
and the factor (k/k')~ w . The latter factor encodes the difference beween Bjorken-x and 
rapidity, y. It is close to one and of minor importance. On the other hand, the factor of k 2 , 
is important as it amplifies the high virtuality contributions to the structure function (note 
that F 2 is obtained from a convolution of the unintegrated gluon density with the photon 
impact factor which selects k 2 values close to the Q 2 of the experiment). 

In Fig. [in]we show the unintegrated gluon density at two values of x - one in the middle 
of HERA region, x = 10~ 3 , and one just above HERA low x region, x = 10 -4 . The 

8 note that the analogy to the QM wave-packet is not complete because in the BFKL eq. rapidity is 
analogous to imaginary time. 
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Figure 13: The unintegrated gluon density using the DP solution of the BFKL equation, 
as determined from the fit to HERA F2 data. 
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Figure 14: Comparison of contributions to the unintegrated gluon density from the eigenfunc- 
tions which are not sensitive (n = 1,2) and which are sensitive to BSM effects (n = 3. .100) 
to the curve at x = 2 x 10~ 3 of Fig. [TT] 
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Figure shows again that the contribution of BSM physics changes the shape of the gluon 
density but does not lead to an increase of the large transverse momentum tail (in fact this 
tail turns out to be smaller for SM+SUSY) contrary to the naive expectation that BSM 
effects should show up as an increased parton activity at large virtualities. The figure also 
shows that the discrepancy between the unintegrated gluon density for the SM alone and 
for SM+SUSY increases substantially with diminishing x at larger k 2 , which suggests that 
Drell-Yan measurements at LHC, in the region of x ~ 10~ 4 and k 2 ~ 50 GeV 2 , could have 
high sensitivity to BSM physics. 

In contrast to the diffusion effect, the important point that we are stressing is that it is 
the spectrum of allowed eigenvalues (and also the shape of the eigenfunctions) that enters 
into the construction of the Green function which is sensitive to physics at high scales. This 
sensitivity emerges from the interplay of two main features of the DP solution: the identity 
of the boundary condition in respect to addition of wk to the phase and the running of the 
coupling constant. For n > 3 this interplay is affected by SUSY effects, whereas for n < 2 
the eigenvalues are unaffected since the critical transverse momentum, t c is below the SUSY 
threshold. For n > 3, the change of the eigenvalues is substantial. As n increases, the 
ratio of eigenvalues with and without SUSY rapidly reaches its asymptotic value of 7/3. In 
Fig. [14] we compare the summed contribution of the eigenfunctions which are sensitive to 
BSM effects (n = 3. ..100) with the sum of contribution of the first two eigenfunctions which 
are not sensitive to SUSY. The contributions of the n = 3. ..100 eigenfunctions changes both 
in size and shape when evaluated with and without SUSY whereas the contribution of the 
first two eigenfunction remains unaltered. 

Fig. [TBI shows the rj variation as a function of the eigenvalue, u, for the fit with the SUSY 
threshold of 10 TeV and for the SM fit only. Both relations show a substantial variation of 
the phase rj with decreasing u. The substantial difference between the two relations reflects 
a large difference between the eigenvalues and eigenfunctions in both cases. 

The changes due to BSM effects lead to a substantial increase of x 2 when the correspond- 
ing gluon densities are confronted with data because the changes are of the order of 10% 
whereas the data precision is around 1 - 2%. Indeed, the evolution of the unintegrated gluon 
density performed with the SM alone - using the same parameter values for t] boundary as in 
the SUSY case (dotted line in Fig. HU dashed lines in Fig. [13]) would give an increase of x 2 
by Ax 2 ~ 160 when compared to data. On the other hand, had we used the r\ boundary pa- 
rameters as determined with the SM alone and performed the evolution with the SM+SUSY 
eigenvalues and eigenfunction, the increase of x 2 would be even larger, A% 2 ~ 300. The dif- 
ferences between the x 2 ' s of the SM and SM+SUSY fits shown in Tabled] are substantially 
smaller, A% 2 ~ 12, because the parameters of the 7] boundary conditions and the proton 
form factors were fitted to data what diminishes substantially the sensitivity to BSM effects. 

It is quite clear that the DP could provide an exciting framework to study BSM effects 
if more can be learned about its infrared boundary. One way to improve our knowledge of 
this boundary is to apply our analysis to an another, independent set of data, such as the 
Drell-Yan processes at the LHC 0. The LHC Drell-Yan data extend in the low-x region to 

9 Natural candidates would be also the diffractive HERA processes. Unfortunately, they are not measured 
with sufficient accuracy because the HERA detectors were not designed to measure diffractive processes. 
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Figure 15: The eigenvalues and infrared phases for the Standard Model and N=l SUSY at 
10 TeV, as determined at k — k . 

much larger Q 2 scales than HERA data. This will allow the study of the evolution in In k 2 
which is dominated in the DP scheme by the low u region. This region and the region of 
very low x is very sensitive to SUSY effects and at the same time is much less sensitive to 
NLO effects which are difficult to compute. 

Future HEP experiments which are now under discussion can also substantially im- 
prove the knowledge of the infrared boundary and of the properties of the DP. The LHeC 
project [28J could provide important information about the region of very low x and not so 
high scales. This could lead to a better understanding of the properties of the large u con- 
tributions. In addition, it will be also possible to measure precisely the exclusive diffractive 
processes. This will provide an independent evaluation of the infrared boundary since in the 
exclusive processes (e.g in the exclusive J /if) or T production) the gluon density contributes 
almost quadratically [2T||23]. 

The properties of the DP can also be studied very well in the 7*7* process in a future linear 
e + e~ collider [29]. The 7*7* process is very interesting because it would permit the direct 
test of the universality of the BFKL pomeron and its boundary conditions. Of particular 
interest could also be the newly proposed electron wakefield accelerator [30], which could 
accelerate electrons into the several TeV energy region. This would allow the measurement 
of the 7*7* structure function at very low x and high Q 2, s, which could further increase the 
sensitivity to BSM effects. 



5 Summary 

In this paper we have analyzed the properties of the Discrete BFKL Pomeron (DP) using 
analytical and numerical methods. We show, using as example N=l SUSY, that BSM physics 
substantially alters its eigenvalue spectrum and the shape of its eigenfunctions. This is a 
genuine sensitivity insofar as it cannot be entirely absorbed into any free parameters of the 
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Discrete Pomeron solution of the BFKL equation. 

The physical origin of this sensitivity can be traced back to the fact that in the low 
Bjorken-x regime the behaviour of the structure functions is directly related to the positions 
of the Regge poles (pomerons). In the BFKL approach, due to the approximate conformal 
invariance, the corresponding eigenfunctions and eigenvalues of the BFKL kernel are deter- 
mined from exponentially large transverse moments, where contributions from particles of 
any BSM physics play an essential role. On the other hand the locality property of the 
BFKL equation assures that in any physical process performed at energy scales which are 
much lower than the BSM ones, the BSM quanta cannot be produced and the transverse 
momenta of the virtual particles are very limited. The effects of large transverse momenta 
appears, however, through the substantial alteration of the eigenstate spectrum of the BFKL 
Hamiltonian, which is important in the low-x region. In our view, this provides a new mech- 
anism for the detection of BSM effects, which has not previously been considered. 

The eigenvalue spectrum of the DP cannot be directly measured because it determines 
the gluon density through a complicated superposition of Pomeron states. The result of 
their interference can be compared with data only after the free parameters of the BFKL 
solution, determining the infrared boundary, are specified. Although the infrared boundary 
is a physical quantity of non-perturbative QCD origin, we could only determine it in this and 
the previous paper within a heuristic approach. Our description of this boundary provides a 
very good fit to the data and shows that, the BSM effects are sensed by the HERA F 2 data, 
notwithstanding the large freedom of the parameter choice. 

The analysis of HERA data indicates an improved quality of fit for the case of N=l 
SUSY, with the SUSY scale as being around 10 TeV. Needless to say that this determination 
is only possible within our heuristic model approach. Our limited knowledge of the infrared 
boundary diminishes substantially (but not completely) the sensitivity of the fit to BSM 
effects. 

This sensitivity can be substantially improved by a better determination of the universal 
boundary condition. We can gain a better understanding of the infrared boundary from 
the analysis of additional data sets, especially of the LHC Drell-Yan data. The data from 
future experimental facilities like LHeC, the e + e~ linear collider or even higher energy plasma 
wakefield accelerators could also become crucial. The 7*7* process which can be very well 
measured at linear colliders is of particular interest since in this reaction the properties of 
the discrete Pomeron solution are simplified owing to the absence of the proton. 

The method described in this paper opens a new possibility of using high precision 
experiments to search for new physics at energy scales considerably larger than the scales 
at which the experiments are performed. We consider the approach formulated here, which 
involves a heuristic model for the parameterization of the infrared phases of the BFKL 
eigenfunctions, as a first step which should be substantially improved by involving more 
data and more theoretical analysis. 
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